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Abstract 

Using an Iterative Leading Particle Model to describing nucleon-air collisions, we determine the hadronic flux 
in the atmosphere and compare with cosmic ray experimental data. 

We present in this paper a calculation of the hadronic flux in the atmosphere. Analytical solutions for 

. . . ^ the nucleonic diffusion equation in the atmosphere having as boundary condition the primary spectrum and 

^ ■ calculated with the leading particle model, shows a strong correlation between the inelastic proton-air cross 

Q>^ . section and the momentum of the leading particle distribution, because the nucleonic shower penetration in 

\ the atmosphere depends on the primary spectrum. In order to analyse experimental data on the hadronic flux 
we correlate the moment of the leading particle distribution in nucleon-air collisions with the respective one 

^ ■ in proton-nucleus collisions using the Glauber model (Glauber 1959; Glauber et al. 1970). 

\ For proton-nucleus scattering, at low energy, several models for describing the leading particle spectrum 

Q ' have been proposed (Iteracting Gluon model and Regge-Mueller formalism) (Yama et al. 1997; Duraes et al. 

CN . 1993; Batista et al. 1998) Here, we shall work in the Iterative Leading Particle Model (Hwa 1984; Hufner & 

I ! Klar 1984) and use the notation of Frichter, Gaisser and Stanev (Fritcher et al. 1997). In this model the leading 

^ \ particle spectrum in p + ^ ^ N(nucleon)+ X collisions is built from sucessive interacions with v interacting 

' proton of the nucleus A and the behaviour is controlled by a straightforward convolution equation. It should 

^ . be mentioned that, strictly speaking, the convolution should be 3-dimensional. Here we only considered the 

lO i 1 -dimension approximation. 

^ ■ It is straightforward to show that in this model the j — th moment of the nucleon-air leading particle 

. distribution, < x'^ >N-air, is correlated wiht the respective nucleon-proton moment, < >N~p, by means 

^ I of the following relation 
p i' 

^.n-"^(l- < >N-a^r) = J dH[l - [^(j) exp[-(l - Ki^)aZT{b)] + 

(l-r?(7)exp[-<r(6)]]] (D 

> ■ 

• 1— I , where T{b) is the nuclear thickness and given by means of the Woods-Saxon model (Woods & Saxon 1954; 
X ■ Barrett & Jackson 1977) 

K„-i(l) = h-i f dyyTSUy) + (1 - fU-i) f dyy ' S;; (y) (3) 

Jo Jo 

The differential nucleonic flux in the atmosphere at some depth t ,is given by 

FNiE,t) = iV„i?-(^+i) exp[-^] (4) 

where NoE~^"'~^^^ is the primary spectrum at t = and A{E) is the nucleonic attenuation length wich is given 
by 

AiE) 24100 ' ^ ^ 
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Figure 1: Nucleonic flux at sea level. Experimental data from (Brook et al. 1964; Ashton et al. 1968; Ashton 
et al. 1970). Continuous line, result of fit. Dash lines, maximal and minimal values of the calculated nucleonic 
flux. 

We tested Eq. Q in comparison with cosmic ray data on nucleonic flux (Brook et al. 1964; Ashton et al. 
1968; Ashton et al. 1970) and hadronic flux in the atmosphere (Mielke et al. 1993; Mielke et al. 1994; Aglietta 
et al. 1997). For a^^f we used the UA4/2 Collaboration parametrization (Burnett et al. 1992) and estimated 
cTj^""*^ by means of the Glauber model (Glauber 1959; Glauber et al. 1970) and for the T{b) nuclear thickness 
we used the Woods-Saxon model (Woods & Saxon 1954; Barrett & Jackson 1977). The parameters r/ and K 
were left free. The best fit /d.o.f. = 2.61) corresponds to r/ = 1 and K = 0.34 and is shown in Fig. 1 
(nucleonic flux at sea level) (Brook et al. 1964; Ashton et al. 1968; Ashton et al. 1970) in Fig. 2 (hadronic flux 
at sea level) (Mielke et al. 1993; Mielke et al. 1994) and in Fig. 3 (hadronic flux at t = SAQg/cw?) (Aghetta 
et al. 1997). The hadronic flux was obtained multiplying the nucleonic flux in Eq. Q by the Kascade factor 
(Mielke et al. 1993; Mielke et al. 1994) 

^^ 7r+ + ^ = 0.04 + 0.27 ln(.E/Gey) (6) 

p+n \ I J 

in order to count the number of pions in the hadronic flux. We have used the same primary spectrum as in 
(Bellandi et al. 1998). We also show in these figures the maximal and minimal values for the flux considering 
the experimental errors in the primary spectrum. 
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Figure 2: Hadronic flux at sea level. Experimental data from (Mielke et al. 1993; Mielke et al. 1994). 
Continuous line, result of fit. Dash lines, maximal and minimal values of the calculated hadronic flux. 




Figure 3: Hadronic flux at i = 830 g/cw?. Experimental data from (Aglietta et al. 1997). Continuous Une, 
result of fit. Dash lines, maximal and minimal values of the calculated hadronic flux. 
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